EVOLUTION OF A UNIVERSE FILLED WITH A CAUSAL VISCOUS FLUID 
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The behaviour of solutions to the Einstein equations with a causal viscous fluid source is in- 
vestigated. In this model we consider a spatially flat Robertson- Walker metric, the bulk viscosity 
coefficient is related to the energy density as £ = ap m , and the relaxation time is given by (/p. 
In the case m — 1/2 we find the exact solutions and we verify whether they satisfy the energy 
conditions. Besides, we study analytically the asymptotic stability of several families of solutions 
for arbitrary m. We find that the qualitative asymptotic behaviour in the far future is not altered 
by relaxation processes, but they change the behaviour in the past, introducing singular instead of 
deflationary evolutions or making the Universe bounce due to the violation of the energy conditions. 



I. INTRODUCTION 



It is believed that quantum effects played a fundamental role in the early Universe. For instance, vacuum polarisation 
and particle production arise from a quantum description of matter. It is known that both of them can be modelled 
OO , in terms of a classical bulk viscosity [l| . 

Cosmological models with a viscous fluid have been studied by several authors. Some of the interesting subjects 
I addressed by them were the effects of viscous stresses on the avoidance of the initial singularity |2J , the dissipation of 
O" 1 a primordial anisotropy Q, the production of entropy Q, and inflation and deflation Q. 

Recently Pavon et al. considered a homogeneous isotropic spatially-flat universe filled with a causal viscous fluid 
bJO whose bulk viscosity is related to the energy density by the power law £ ~ p m [6] . 

In this paper we extend and improve the analysis of this cosmological model. We present the set of equations 
that describe the model in section 2. In section 3 we solve them when to = 1/2 and give a detailed analysis of their 
solutions. In section 4 we study of stability of those asymptotic solutions that appear for m ^ 1/2 in the noncausal 
OO ' model by means of the Lyapunov method. The conclusions are stated in section 5. 

m 

II. THE MODEL 

OO 

© : In the case of the homogeneous, isotropic, spatially flat Robertson-Walker metric 



ds 2 = dt 2 - a 2 {t){dx\ + dx 2 2 + dxl,) (1) 

only the bulk viscosity needs to be considered. Thus we replace in the Einstein equations the equilibrium pressure p 
by an effective pressure [3] 

H 2 = \p H + 3H 2 = l{p- p - a ) (2) 

where H= a/a, ' — d/dt, p is the energy density, a is the viscous pressure, and we use units c = 8irG = 1. As equation 
of state we take 

p = ( 7 -l)p (3) 
with a constant adiabatic index 7 > 0, and a has the constitutive equation of a viscoelastic fluid 

0- + TO- = -3(H (4) 

Here £ > is the bulk viscosity coefficient and r is the bulk relaxation time and causality demands r > 0. Following 
we choose 
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C = ap m r = C/p (5) 
where a and m are constants. Using Eqs. 2-5 we get 

3 HI \ ' HI ) 2' \ HI j ' 7- , it- K J 

^ J ff + (e + 77o)^+^(£7-7o)ff 3 = ! r = (7) 

H+*Lll\H\ 1 -*-lH a = > r^O, 7 = (8) 
where e = sgnH, r = 2m — 1, H = (j/jo) 1 ^ /V3 and 70 = \fZa. We assume a > 0. 

III. CASE r = 

The general solution of Eq. [7] for H > takes the following parametric form: 

l + 77o 1 + 77o J « W) 

where A± are the roots of A 2 + A + J°^ a °l = and A, B are arbitrary integration constants. 

A. One-Parameter Solutions 

The one-parameter families of solutions arise when either A or B vanishes and can be obtained explicitly: 

H±(t) = i/±/At, 7^70 (10) 
H+ = D, H_=u /At, 7 = 70 (11) 



v = 1 + 770 , ! 
± 3(7-70) 



_ 470(7-70) 
(1 + 77o) 2 



" = 3rT^ (12) 



Thus, expanding solutions '-' are always Friedmann, but solutions '+' are Friedmann for 7 > 7o , de Sitter for 7 = 70 
and explosive for 7 < 70. 

B. Two-Parameter Families of Solutions 

When AH ^ 0, a(t) can be written in closed form in terms of known functions only for some values of 7 and 
70. However, in general, we need to study the solution in the parametric form of Eq. [9] . We obtain the following 
classification for the two-parameter families of solutions : 

A. The evolution occurs between singularities, it reaches a maximum and recollapses again. The leading behaviour 
near the singularities is Friedmann, as 1101 - for 7 ^ 70 or 1111 - if 7 = 70 ■ These solutions have particle horizons 
(0 < is- < 2/3). 

B. There is a bounce with an explosive singularity in the past and its behaviour in the future is either: 

1. asymptotically Friedmann, as HON - for 7 > 7 o- 

2. asymptotically de Sitter, asHTlf for 7 = 70. 

3. divergent at finite times, with leading behaviour [TOlf for 7 < 70. 

C. The evolution begins at a singularity with a Friedmann leading behaviour, as 1101 - for 7 7^ 70 or 1111 - if 7 = 70; and 
so they have also particle horizons. Then it expands, and its behaviour in the future is like B. 

D. The evolution begins at an explosive singularity and ends at a big-crunch singularity. 



3 



C. Energy Conditions 

• Dominant energy condition (DEC): p > \p + a\ —3H 2 < H < 0. 

DEC is violated part of the time in the following families: A about the maximum (if v_ > 1/3); Bl, about the bounce 
as well as for large times (if v + > 1/3); near the singularity in CI, C2 (if v < 1/3 ) and C3 (if v_ < 1/3); C3, near 
the "explosion". DEC is violated always in families A (if < 1/3), B2 and B3. 

• Strong energy condition (SEC): p + 3p + 3a > H + H 2 < 0. 

SEC is satisfied always in families A and CI (if v + < 1). It is violated part of the time in the families: Bl, about the 
bounce (if i/+ < 1); for large times in CI (if v + > 1), C2 and C3. SEC is violated always in families Bl (if v + > 1), 
B2 and B3. 



IV. CASE r # 



We investigate the asymptotical stability of behaviors that occur in the noncausal model by means of the Lyapunov 
method. 



A. Stability of the de Sitter solution 

For the de Sitter solution H = Hq we rewrite Eq. [S] as a " mechanical system" 
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H + V(H) 
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= -3H 2 H 


7 + - 
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(13) 



Q Q HZ 

V(H) = --H 4 + --^H^ r ± 4 
v ' 8 24-r T 



(14) 
(15) 



B. Stability of the Asymptotically Friedmann solution 

For r > it is easy to check that Eq. [S] admits a solution whose leading term is 2/(37*). To study its stability we 
make the change of variables H = v(z)/t, t r — z. Then this equation takes the form 



d 

dz 



l -v' 2 + V(v,z) 



(16) 



4 — r 37(3 — r) 



v 3 ~ r „ ( 1 
z \ z z 



(17) 



V{v, Z ) = \Hl(v-l-lnv)\ + o(±). r = A 



(18) 



Stability of exponential-like solutions 



In the case of the solution for 7 = and m 7^ 

2m 



i(t) = exp 



2m- 1 



w M t~ 



W M 



= (3 ro+1 ma)" 



(19) 
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we make the change of variables H = (sw(s)/t), t r /( r + 1 ) = s and get 



d_ 

d,s 



^w' 2 + W(w, s) 



where now the potential to leading order in 1/s is 



3 l-r/2 r+1 w ,2 

-s + Ol - 



7o r w 



r-l ' 



(20) 



r + 1 



3 l-r/2 w 3-r 0,( r + l) z 

(3-r) 70 + ~ 



r + 3 



(21) 



gl-r/2 

s) = — 4— lnu; — u; 4 , r = 3 

97o 



(22) 



CONCLUSIONS 



When r = 0, the splitting of the large time asymptotic behavior of solutions in terms of sgn (7 — 70 ), closely 
resembles the classification for the noncausal solutions. However, causality makes new families of solutions appear as 
the bouncing ones and those which expand from a singularity but recollapse in a finite time at another singularity. 
Most singular solutions have particle horizons. 

We demonstrate that there is no value of r for which there is a stable expanding de Sitter period in the far past. 
This supports strongly the conclusion of that causality avoids the deflationary behavior proposed by Barrow Q . 

If r < 0, a stable inflationary phase occurs in the far future for any 7 > 0; the condition 7 = 70 is required if r = 0, 
and such a behavior is unstable for r > 0. We note that the noncausal model has no stable de Sitter solution if r = 0. 
If 7 = and m < 0, the same faster than exponential expansion found by Barrow Q is asymptotically stable. 

If r > 0, we find that relaxation effects do not alter the perfect fluid behavior a ~ i 2 /( 37 ) for t — > 00. This arises 
because the viscous pressure decays faster than the thermodynamical pressure. However, if r = and 7 > 70, both 
pressures decay asymptotically as t~ 2 and the exponent becomes v + . The perfect fluid behavior becomes unstable if 
r < 0. 

Large negative viscous pressures may arise in stable evolutions which avoid an initial singularity or have a viscosity- 
driven inflationary stage. If r = 0, no two-parameter solution satisfies the energy conditions. 
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